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A  dynamic  model  of  a  discrete  reversible  fuel  cell  (RFC)  system  has  been  developed  in  a  Matlab  Simulink® 
environment.  The  model  incorporates  first  principles  dynamic  component  models  of  a  proton  exchange 
membrane  (PEM)  fuel  cell,  a  PEM  electrolyzer,  a  metal  hydride  hydrogen  storage  tank,  and  a  cooling  system 
radiator,  as  well  as  empirical  models  of  balance  of  plant  components.  Dynamic  simulations  show  unique 
charging  and  discharging  control  issues  and  highlight  factors  contributing  to  overall  system  efficiency. 
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1.  Background  and  introduction 

Reversible  fuel  cell  (RFC)  systems  act  as  “rechargeable”  energy 
storage  devices  designed  to  work  as  turn-key  systems  with  electric¬ 
ity  being  the  only  input  and  output.  This  may  allow  future  systems 
to  replace  traditional  chemical  storage  batteries  in  applications 
where  the  potential  performance  benefits  outweigh  the  cost  trade¬ 
offs  [1  ].  The  major  components  in  an  RFC  system  are  a  fuel  cell,  an 
electrolyzer  ( unitized  RFC  systems  use  just  one  device  to  perform 
both  functions,  discrete  systems  have  separate  fuel  cell  and  elec¬ 
trolyzer  devices),  and  a  hydrogen  storage  tank.  Numerous  balance 
of  plant  components  are  also  necessary  for  a  practical  system  [2]. 
It  is  only  through  thorough  modeling  of  all  these  components  that 
true  system  efficiencies  can  be  estimated.  Fig.  1  shows  a  schematic 
diagram  of  the  RFC  system  modeled  in  this  work,  demonstrating 
the  complex  auxiliary  components  and  circuitous  integration  of  all 
of  the  parts. 

The  RFC  system  is  “charged”  by  supplying  electricity  to  an  elec¬ 
trolyzer  that  disassociates  water  yielding  hydrogen  and  oxygen.  The 
oxygen  is  vented  to  the  atmosphere  and  the  hydrogen  is  stored 
in  the  hydrogen  tank.  The  current  system  uses  a  reversible  metal 
hydride  compound  confined  in  an  aluminum  tank  for  hydrogen 
storage.  When  power  is  needed  from  the  system,  hydrogen  is  sup¬ 
plied  by  the  tank  to  the  fuel  cell  where  it  is  combined  with  air  to 
produce  electricity,  heat,  and  water. 
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The  RFC  energy  storage  system  is  unique  when  compared  to  tra¬ 
ditional  chemical  batteries  because  the  power  capability,  energy 
storage  capacity,  and  recharge  rate  are  all  determined  indepen¬ 
dently  [3].  With  the  exception  of  shared  controls  and  plumbing, 
the  hydrogen  tank  capacity  solely  determines  the  energy  storage  of 
the  system,  the  fuel  cell  power  output  governs  the  power  output  of 
the  system,  and  the  electrolyzer  power  determines  the  “recharge” 
rate  of  the  system.  These  properties  lead  to  system  advantages 
in  specialized  applications,  particularly  for  those  applications  that 
require  low  to  moderate  power,  long  duration  performance,  and/or 
strict  weight  requirements.  The  fuel  cell  and  electrolyzer  can  be 
sized  as  needed  for  the  power  requirement,  while  the  hydrogen 
storage  tank  can  simply  be  sized  up  or  down,  independently,  to 
meet  the  energy  storage  requirement.  A  traditional  chemical  bat¬ 
tery  would  require  increases  in  all  parameters  in  order  to  meet  the 
demands  of  just  one  of  these  design  criteria.  Consequently,  a  fuel 
cell/electrolyzer  system  can  theoretically  achieve  a  better  energy 
density  than  even  state-of-the-art  chemical  batteries  [4]. 

In  order  to  evaluate  the  performance  of  a  fuel  cell  rechargeable 
energy  system,  to  garner  insight  into  the  dynamic  response  char¬ 
acteristics  of  the  system  and  individual  components,  and  to  design 
superior  systems  in  the  future,  a  dynamic  model  of  a  RFC  system  is 
developed  in  a  modular  Matlab  Simulink®  framework. 

t.l  Experimental  RFC  system  description 

The  particular  regenerative  fuel  cell  system  being  modeled  was 
designed  to  replace  lead-acid  batteries  onboard  a  U.S.  Army  Stryker 
vehicle  and  consists  of  a  5.5  kW  fuel  cell,  a  3.6  kW  electrolyzer,  and 
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Fig.  1.  Schematic  representation  of  the  experimental  RFC  system. 


a  hydride  tank  capable  of  holding  10,000  standard  liters  of  hydro¬ 
gen.  This  combination  of  components  is  assembled  using  parts  that 
are  commercially  available,  with  little  effort  made  to  optimize  the 
system  for  any  particular  performance  goal.  As  a  result,  the  perfor¬ 
mance  reported  herein  should  not  be  construed  as  optimal.  Future 
system  performance  improvements  are  very  likely  through  better 
component  sizing,  integration,  and  optimization,  which  is  outside 
the  scope  of  this  paper. 

The  Stryker  vehicle  is  equipped  with  a  large  array  of  advanced 
electronic  surveillance  equipment,  which  can  be  utilized  in  a 
“silent-watch”  mode  whereby  the  electronics  are  powered  for  an 
extended  time  without  running  the  diesel  prime  mover.  The  system 
is  charged  by  excess  power  from  the  vehicle’s  alternator  dur¬ 
ing  normal  operation,  but  when  silent  operation  is  desired,  the 
main  vehicle  combustion  engine  is  turned  off,  the  vehicle  becomes 
stationary,  and  the  batteries  power  the  onboard  electronics.  The 
proposed  RFC  system  would  seamlessly  replace  the  lead-acid  bat¬ 
tery  pack  while  ideally  providing  longer  silent-watch  duration  due 
to  the  efficiency  and  energy  storage  capabilities  of  the  RFC  system. 

RFCs  are  well  suited  to  military  applications  because  of  the 
desire  to  use  common  fuels.  Jet  Propulsion  8,  or  JP-8,  is  a  kerosene- 
based  fuel  that  is  the  single  battlefield  fuel  for  Department  of  the 
Army  applications  including  electric  power  generators,  wheeled 
and  tracked  vehicles,  aircraft,  stoves,  and  heaters.  An  immense 
logistical  effort  is  necessary  to  move  the  enormous  quantities  of 
JP-8  fuel  required  for  military  operations  into  a  forward  deployed 
overseas  location.  Delivered  costs  of  JP-8  to  Army  combat  platforms 
have  been  estimated  to  be  at  least  $30-$40  gallon-1  for  overland 
transport,  and  greater  than  $400  gallon-1  for  air  delivery  [5]. 

Because  of  the  “single  battlefield  fuel”  mandate,  fuel  cell  gen¬ 
erators  or  other  clean  and  efficient  hydrogen  consuming  devices 
cannot  be  utilized  on  the  battlefield  unless  the  hydrogen  is  obtained 
from  reformation  of  JP-8  or  by  other  means  with  existing  battle¬ 
field  resources.  Consequently,  an  RFC  system  having  electricity  as 
the  only  input  and  output  is  one  potential  fuel  cell  technology  avail¬ 
able  to  the  military.  This  type  of  RFC  device  could  be  recharged  by 
alternators  attached  to  diesel  engines  in  exactly  the  same  manner 
batteries  are  used  today.  This  could  benefit  military  applications 
that  currently  rely  on  lead-acid  batteries  for  electrical  storage,  with¬ 
out  violating  the  single  fuel  forward  mandate. 

A  proton  exchange  membrane  (PEM)  fuel  cell  and  electrolyzer 
were  chosen  for  this  application  because  of  their  relatively  high 
power  density,  inherent  load  following  capability,  and  relative  tech¬ 
nological  maturity  compared  to  other  fuel  cell  and  electrolyzer 


technologies.  The  fuel  cell  is  fueled  by  pure  hydrogen  produced 
from  electrolysis  of  deionized  water  and  stored  in  a  metal  hydride 
tank.  Metal  hydride  absorbs  and  desorbs  hydrogen  in  exothermic 
and  endothermic  reactions,  respectively  [6].  Consequently,  the  RFC 
system  utilizes  an  integrated,  dual-purpose,  cooling  loop  that  pro¬ 
vides  coolant  to  the  hydride  tank  during  the  exothermic  charging 
process,  and  fluid  warmed  by  the  fuel  cell  waste  heat  to  the  tank 
during  the  endothermic  hydrogen  release  process.  The  use  of  a 
metal  hydride  tank  therefore  creates  a  synergy  that  enables  fuel  cell 
waste  heat  to  provide  energy  for  hydrogen  storage,  eliminating  the 
need  for  compression  energy  that  would  be  required  if  compressed 
gas  storage  was  used. 

12.  RFC  models  in  the  literature 

Several  RFC  system  models  and  analyses  have  been  reported 
in  the  literature.  Bolwin  calculated  static  RFC  system  energy  and 
power  densities  based  on  a  list  of  component  masses,  efficien¬ 
cies,  and  power  densities  [3].  Ulleberg  developed  a  semi-empirical, 
bulk,  alkaline  electrolyzer  model  for  use  in  hydrogen  energy  system 
simulations  and  showed  an  example  system  integration  including 
renewable  energy  sources  and  compressed  gas  hydrogen  storage 
[7].  Barbir  et  al.  statically  compared  the  efficiencies  of  unitized  and 
discrete  PEM  RFCs  using  generalized  polarization  curves  and  par¬ 
asitic  losses  proportional  to  system  power  [4].  Busquet  et  al.  used 
an  empirical  model  to  well  predict  fuel  cell  and  electrolyzer  per¬ 
formance,  but  did  not  include  parasitic  ancillary  loses  or  energy 
storage  [8].  Khan  and  Iqbal  created  a  dynamic  Simulink  model  of 
a  wind  powered  fuel  cell  electrolyzer  system,  but  did  not  model 
parasitic  loads  or  heat  transfer  within  the  devices  [9].  Gorgiin  pro¬ 
duced  a  detailed  dynamic  PEM  electrolyzer  model  for  future  use  in 
fuel  cell  electrolyzer  systems  that  includes  pressurized  hydrogen 
storage  and  membrane  water  transport.  Flowever,  the  model  does 
not  include  heat  transfer  or  GDL  gas  diffusion  effects  [10].  Maclay  et 
al.  developed  a  comprehensive  empirical  dynamic  residential  RFC 
model  used  to  analyze  RFC  applications  to  renewable  energy  stor¬ 
age,  though  the  analysis  does  not  include  efficiency  losses  due  to 
ancillary  components  or  hydrogen  storage  effects  [11].  The  model 
reported  herein  differs  from  others  in  the  literature  due  to  the 
exhaustive,  first  principles  approach  that  has  been  undertaken.  This 
approach  enables  simulation  of  the  dynamic  responses  of  all  the 
major  system  components  and  all  of  the  interactions  amongst  the 
components  and  in  the  balance  of  plant,  along  with  comprehensive 
energy  and  efficiency  analyses  that  include  all  significant  ancil- 
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lary  components  of  a  complete  RFC  system.  The  model  is  detailed 
enough  to  accurately  simulate  this  type  of  system,  yet  simple 
enough  to  use  in  control  system  development  and  analysis  studies. 

2.  Model  description 

The  RFC  model  simulates  the  experimental  RFC  system 
described  above  as  closely  as  possible  by  requiring  only  the  same 
inputs  as  provided  to  the  experimental  system.  Operating  input 
parameters  include  ambient  air  temperature  and  humidity,  desired 
operating  mode  (power  mode  or  charge  mode)  and  desired  fuel 
cell  current  when  in  power  mode.  Extensive  material  property 
data,  component  geometry,  and  dimensional  information  are  also 
entered  into  the  model  to  simulate  the  properties  of  the  experi¬ 
mental  system  components. 

The  model  assumes  that  maximum  power  is  always  available 
to  the  electrolyzer  when  it  is  operating  in  charge  mode.  This 
simplifying  assumption  allows  the  dynamics  of  the  system  to  be 
analyzed  separately  from  the  electrical  power  source  (an  alterna¬ 
tor  attached  to  a  diesel  engine  in  the  current  experiment).  This 
assumption  would  not  be  valid  in  all  circumstances;  for  example, 
if  the  electrolyzer  was  powered  by  an  intermittent  source  such  as 
a  photovoltaic  panel,  the  maximum  power  would  not  always  be 
available.  Such  a  case  would  yield  interesting  results,  but  may  give 
little  information  about  the  performance  of  the  RFC  system  itself, 
which  is  of  primary  interest  in  the  current  effort. 

Despite  the  attempt  to  make  this  model  as  comprehensive  as 
possible,  numerous  balance  of  plant  components  that  are  neces¬ 
sary  for  a  real  system  to  work  have  been  assumed  to  have  no  impact 
on  the  dynamic  performance  of  the  system.  These  balance  of  plant 
components  that  are  assumed  to  have  instantaneous  response  char¬ 
acteristics  mainly  include  plumbing  items  such  as  valves,  piping, 
and  filters.  For  example,  the  deionized  water  system  necessary  for 
robust  electrolysis  has  been  neglected  and  flow  pressure  losses  have 
not  been  accounted  for  in  any  part  of  the  model.  In  addition,  no 
effort  has  been  made  to  produce  a  physical  model  of  the  enthalpy 
and  humidity  exchange  between  fuel  cell  cathode  outlet  and  the 
fuel  cell  inlet.  While  these  omissions  affect  system  performance  the 
major  contributors  to  overall  system  performance  are  captured.  The 
effects  of  these  missing  balance  of  plant  parts  could  be  introduced 
into  the  model  in  a  later  study,  but  are  assumed  to  be  minimal  in 
the  current  analyses. 

Each  of  the  main  components  of  the  RFC  model  (fuel  cell,  elec¬ 
trolyzer,  metal  hydride  tank,  radiator)  is  simulated  using  a  dynamic, 
first  principles  approach.  The  environments  within  any  control 
volume  of  the  models  are  assumed  to  comprise  ideal,  “perfectly 
stirred”  conditions  for  fluids  and  homogenous  compositions  for 
solids.  The  following  sections  will  describe  the  conservation  equa¬ 
tions  and  physics  modeled  for  each  of  the  four  main  components. 
The  equations  and  approach  used  for  each  of  the  models  are  sim¬ 
ilar,  but  the  actual  terms  and  physical  properties,  and  processes 
simulated,  vary  with  each  component  model.  Table  1  shows  the 


Table  1 

Control  volumes  and  physical  relations  modeled  for  each  of  the  four  major  system 
components 


Model  attribute 

Fuel  cell 

Electrolyzer 

Hydride  tank 

Radiator 

Number  of  control  volumes 

8 

6 

31 

8 

Conservation  of  energy 

X 

X 

X 

X 

Conservation  of  mass 

X 

X 

X 

X 

Heat  conduction 

X 

X 

X 

X 

Heat  convection 

X 

X 

X 

X 

Diffusion 

X 

X 

Chemical  reactions 

X 

Electrochemical  reactions 

X 

X 

number  of  control  volumes  and  the  particular  physics  modeled  for 
each  component. 

2  A.  Metal  hydride  tank  model 

The  metal  hydride  tank  is  modeled  as  a  cylinder  comprised 
of  30  nested  cylindrical  shells  as  shown  in  Fig.  2.  The  hydrogen 
gas  can  enter  or  exit  from  only  one  end  of  the  tank.  In  addition 
to  the  reactive  material  “shells”  pictured  in  Fig.  2,  a  thin  alu¬ 
minum  outer  enclosure  is  also  modeled  that  does  not  participate 
in  the  hydriding  reaction,  but  simulates  the  canister  that  holds 
the  powdered  material  in  an  actual  tank.  Note  that  heat  transfer 
physics  are  affected  by  the  thermal  mass,  conduction,  and  con¬ 
vection  heat  transfer  properties  of  the  aluminum  tank.  Previous 
work  [12]  has  shown  that  accurate  tank  temperature  and  hydrogen 
absorption/desorption  results  can  be  obtained  from  a  metal  hydride 
model  without  accounting  for  gaseous  mass  transport  within  the 
hydride  bed  because  the  timescales  of  hydrogen  diffusion  through 
the  porous  hydride  material  are  several  orders  of  magnitude  faster 
than  the  physics  associated  with  heat  transfer  within  the  tank.  As 
a  result,  hydrogen  absorption  and  desorption  are  rate-controlled 
by  heat  transfer  physics  [12].  Likewise,  it  has  been  shown  that 
gaseous  heat  transfer  and  gaseous  energy  conservation  can  also  be 
neglected  as  the  quantity  of  energy  transferred  to  or  from  the  gas  is 
several  orders  of  magnitude  smaller  than  the  energy  of  the  hydrid¬ 
ing  reaction,  on  a  per  mass  basis  [12].  These  simplifications  enable 
the  relatively  large  number  of  30  nodes  to  be  modeled  without 
creating  an  overwhelming  computational  burden.  A  minimum  of 
30  nodes  was  found  to  be  necessary  to  accurately  predict  the  heat 


Not  drawn  to  scale 


Fig.  2.  Drawing  showing  nodal  discretization  used  in  the  metal  hydride  tank  model. 
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transfer,  and  consequently  the  hydrogen  absorption  level,  in  the 
tank  due  to  the  heavy  dependence  of  absorption  on  temperature 
and  the  relatively  poor  heat  transfer  characteristics  of  powdered 
metal  hydride. 

No  absorption-desorption  hysteresis  is  accounted  for  in  the 
present  model.  The  hysteresis  associated  with  the  LaNi5  alloy  used 
herein  is  0.13,  meaning  that  13.8%  more  pressure  is  required  to  fill 
the  tank  than  is  returned  when  it  is  emptied  [6].  This  effect  would 
serve  to  lower  the  overall  RFC  system  efficiency. 

2.1.1.  Rate  of  hydride  reaction 

The  hydriding  reaction  rate  is  modeled  as 

^reaction  —  — ^a  ^  ^  In  (Pfull  —  Aiode)^S  (1) 

where  Ca  is  the  reaction  rate  constant,  Ea  is  the  activation  energy, 
R  is  the  ideal  gas  constant,  and  pfun  and  pno de  are  the  saturated 
and  nodal  densities  of  the  solid  alloy,  respectively  [13].  Pg  is  the 
hydrogen  pressure  within  the  node,  and  due  to  the  fast  transport 
of  hydrogen  through  the  powered  metal,  Pg  is  always  assumed  to 
be  equal  to  the  head  pressure  at  the  inlet  to  the  tank.  Pvh  is  the 
pressure  predicted  by  the  van’t  Hoff  equation  [6]: 

,  ,  AH  AS 

11  vH  ~~  RTg  fT  ^ 

The  values  of  enthalpy  of  reaction,  AH,  and  the  change  in  entropy 
for  the  reaction,  A5,  were  experimentally  determined  for  the  given 
hydride  alloy  being  used  [6]. 

2.1.2.  Solid  alloy  conservation  of  mass 

The  volume  of  solid  hydride  alloy  within  a  node,  Vs.  is  equal  to 
the  non-porous  volume  of  the  node: 

VS  =  (l-£)Vnode  (3) 

where  £  is  the  porosity  of  the  metal  hydride.  The  mass  of  the  solid, 
ms,  in  the  node  is 

—  Pnode^S  (4) 

The  mass  in  a  node  can  only  change  by  absorbing  or  desorbing 
hydrogen.  The  mass  conservation  equation  is  consequently  written 
as 

dms  _  . 

^  —  ^reaction  (pj 

This  can  be  integrated  to  find  the  instantaneous  nodal  mass.  When 
plugged  back  into  Eq.  (4),  this  gives  the  value  of  metal  hydride  den¬ 
sity  needed  to  calculate  the  rate  of  hydrogen  absorption  given  by 
Eq.(l). 

2.1.3.  Conservation  of  energy 

The  change  in  energy  of  the  solid  alloy  in  each  node  is  due  to 
three  processes:  (1)  conduction  with  solid  in  surrounding  nodes, 

(2)  heat  release  or  generation  during  the  sorption  reaction,  and 

(3)  mass  change  via  hydrogen  sorption.  The  energy  conservation 
equation  for  the  solid  within  each  node  is  therefore: 

dEs 

—  Qconduction  +  ^reaction  ^H  (6) 

where  the  rate  of  change  of  solid  energy  can  be  written  as 

dEs  _  ^(PnodeVsChTs) 

dt  ~  at  {  } 

The  nodal  solid  temperature  can  be  found  by  combining  (6)  and  (7) 
and  numerically  integrating. 

As  mentioned  previously,  for  computational  efficiency  and  neg¬ 
ligible  error,  the  temperature  of  the  hydrogen  gas  in  each  node 


Table  2 

Metal  hydride  tank  parameters 


Parameter 

Value 

Units 

Reference 

Number  of  tanks 

3 

Per  design 

Hydride  material 

LaNi5 

Per  design 

Canister  material 

Aluminum 

Per  design 

Canister  thickness 

1.6 

mm 

Per  design 

Total  hydrogen  capacity 

10,000 

SL 

Per  design 

Maximum  pressure,  Pmax 

689.5 

kPa 

Per  design 

Length,  l 

0.384 

m 

Per  design 

Diameter,  d 

0.148 

m 

Per  design 

Hydride  conduction 
coefficient,  kh 

1.0 

W(mK)-1 

[14] 

Aluminum  conduction 
coefficient,  kM 

237 

W(mK)-1 

[15] 

Hydride  specific  heat,  Ch 

418.7 

J(kgK)-1 

[16] 

Al  specific  heat,  CAi 

903 

J(kgK)-1 

[15] 

Density  of  Al,  pAi 

2700 

kg  m-3 

[15] 

Density  of  hydride  metal, 

Ph 

8300 

kg  m-3 

[6] 

Solid  porosity,  e 

0.44 

Calculated 

Activation  energy,  £a 

31,000 

J  mol-1 

[14] 

Reaction  rate  constant,  Ca 

2800 

s-1 

[17] 

Enthalpy  of  reaction,  AH 

30,800 

J  mol-1 

[6] 

Entropy  of  reaction,  AS 

108 

J(molK)-1 

[6] 

Convection  coefficient 
between  canister  and 
coolant,  h 

700 

W(m2  K)-1 

[12]a 

a  This  value  was  established  experimentally  for  a  cylindrical  tank  subjected  to  a 
circulating  water  bath. 


is  assumed  to  be  the  same  as  the  temperature  of  the  solid  metal 
hydride  in  that  node.  Consequently,  in  all  parts  of  the  tank: 

Tg  =  Ts  (8) 

2.14.  Heat  transfer 

Heat  transfer  within  the  tank  is  governed  by  Fourier’s  law: 
Qconduction  —  (9) 

where  A  is  the  cross-sectional  area  between  adjacent  nodes  and  x 
is  the  distance  between  nodal  centers.  Eq.  (9)  applies  to  conduction 
between  the  metal  hydride  nodes  and  between  the  outermost  node 
and  the  aluminum  enclosure.  Convective  heat  transfer  described  by 
Newton’s  Law  of  Cooling  takes  place  between  the  outer  surface  of 
the  aluminum  enclosure  and  the  circulating  coolant  as  follows: 

Qconvection  —  EAC  d T  (10) 

Here  Ac  refers  to  the  outer  surface  area  of  the  aluminum  canister. 
2.1.5.  Model  parameters 

The  system  modeled  in  this  paper  uses  three  individual,  cylin¬ 
drical  tanks  cooled  by  submersion  in  a  flowing  water  bath,  which 
approximates  the  experimental  configuration.  The  specifications 
for  these  tanks  are  given  in  Table  2. 

2.2.  PEM  fuel  cell  model 

The  PEM  fuel  cell  model  used  in  this  system  model  was  devel¬ 
oped  by  Mueller  et  al.  at  UC  Irvine  and  is  described  completely  in  a 
multi-nodal  form  in  [  18  ].  The  bulk  version  of  the  model  used  herein 
is  similar,  with  the  primary  differences  being  the  number  of  nodes, 
the  number  of  cells  in  the  stack,  and  the  dimensions  of  the  different 
parts  of  the  stack.  The  current  first  principles  PEM  fuel  cell  model 
resolves  conservation  of  mass,  conservation  of  energy,  and  species 
conservation  for  each  species  and  layer  within  a  repeat  cell  unit  of 
the  stack. 
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Fig.  3.  Diagram  showing  eight  control  volumes  modeled  in  the  PEM  fuel  cell  com¬ 
ponent. 


The  model  is  divided  into  eight  control  volumes:  anode  and 
cathode  gas  diffusion  layers,  bulk  gas  chambers,  current  collectors, 
membrane,  and  cathode  side  coolant  channel,  as  shown  in  Fig.  3. 
The  first  principles  approach  enables  prediction  of  species  concen¬ 
trations  and  temperatures  for  each  gas  or  solid  in  each  layer  of  the 
cell  repeat  unit  and  the  output  cell  voltage  at  any  time.  The  stack 
outputs  are  simply  determined  by  multiplying  or  adding  the  cell 
repeat  unit  outputs  as  determined  by  the  number  of  series  or  par¬ 
allel  cells  in  the  stack.  For  example,  for  unit  cells  connected  in  series 
the  total  stack  voltage  is  given  by 

^total  —  ft^cell  (11 ) 

where  n  is  the  number  of  cells  in  series.  In  this  case  the  total  stack 
current  is  the  unit  cell  current.  Similar  calculations  are  used  to 
determine  total  stack  species  flux  values. 

The  multi-nodal  model  reported  in  [18]  has  been  validated 
against  actual  operating  data  from  several  PEM  fuel  cells.  How¬ 
ever,  due  to  the  large  computational  requirements  of  the  complex 
multi-nodal  model,  the  reduction  to  just  one  node  has  been  made 
to  allow  integration  into  the  larger  RFC  system  model.  Experience 
shows  that  this  model  still  produces  good  predictions  of  stack  tem¬ 
perature,  voltage,  and  outlet  species  concentrations. 


2.2.1.  Mass  conservation 

Mass  must  be  conserved  in  the  bulk  gas,  GDL,  and  membrane 
sections  of  the  fuel  cell  model  while  tabulating  the  proper  species 
concentrations  of  hydrogen,  oxygen,  nitrogen,  and  water.  The  gen¬ 
eral  species  conservation  equation  for  a  control  volume  can  be 
written  as 

^  =  iiin  -  iiout  +  lireact  +  ^  "diffuse  (general  form)  (12) 

where  hin  and  hou t  are  the  convective  molar  fluxes  in  and  out  of 
the  control  volume,  hreact  is  the  molar  rate  of  species  production, 
and  ^2  hdiffuse  is  the  sum  of  the  species  diffusion  rates.  This  is  the 
most  general  mass  conservation  equation.  Certain  terms  drop  out 
depending  upon  the  particular  control  volume  that  is  being  mod¬ 
eled.  For  example,  only  the  reaction  and  diffusion  terms  are  used 
in  the  GDL  layers  where  gases  diffuse  and  react: 

^  =  "react +  ^  "diffuse  (GDL)  (13) 


Only  diffusion  and  bulk  convective  flow  terms  are  used  in  the  bulk 
gas  layers  because  no  reactions  are  occurring: 

^  =  "in  -  " out  +  "diffuse  (bulk  gas  layers)  (14) 


And  the  membrane  species  conservation  equation  contains  only 
diffusion  terms  with  no  bulk  fluid  flow  or  reaction: 

^  ^  "diffuse  (membrane)  (15) 


Species  conservation  is  maintained  by  integrating  these  equations 
for  each  species  in  each  control  volume  within  the  model. 

For  each  control  volume,  relative  humidity  is  calculated  as  a 
function  of  temperature.  Any  mass  of  water  exceeding  saturation  is 
modeled  as  condensed  liquid. 


2.2.2.  Energy  conservation 

The  general  conservation  of  energy  equation  used  in  each  fuel 
cell  control  volume  is  written  as 


5>inhin 

—  ^  >  thouthout  E  Qht  +  ^  Qreact  (§eneral  form) 

(16) 


where,  on  the  left  side,  m  is  the  mass  of  the  component,  Cv  is  the 
constant  volume  specific  heat,  T  is  temperature,  and  t  is  time.  On 
the  right,  rhin  and  fnout  are  the  mass  transfer  rates  entering  and  exit¬ 
ing  the  control  volume  and  hin  and  hou t  are  the  associated  enthalpy 
values,  respectively.  Heat  transfer  with  surrounding  nodes  is  given 
by  Qht  and  any  heat  generated  or  consumed  by  chemical  reac¬ 
tions  is  given  by  Qre act-  Again,  the  terms  used  depend  upon  the  fuel 
cell  component  being  analyzed.  For  the  current  collectors,  Eq.  (16) 
reduces  to 

msoiidC^- =  ^^Qht  (current  collectors)  (17) 

because  the  solid  components  have  no  enthalpy  transferred  by  con¬ 
vection  nor  any  reactions  occurring  within  the  solid.  For  the  MEA 
(both  GDL  layers  and  the  membrane  are  assumed  to  be  at  the  same 
temperature)  all  of  the  terms  in  Eq.  (16)  are  used  with  the  caveat 
that  the  enthalpy  flow  term  refers  to  energy  transferred  via  diffu¬ 
sion,  not  bulk  convective  flow: 

f  ^  >  tTlgasCy  +  ^  ^  TtlsolidC^)  ^  ^  ^diffuse  ^diffuse  T  ^  ^  Qht 

+  ^&eact  (MEA)  (18) 

For  the  bulk  gas  layers,  energy  is  transferred  by  heat,  bulk  convec¬ 
tive  flow,  and  diffusion,  but  no  reaction  energy  is  included: 

^  HlgasCv  +  ^  ^  ^solid^  ^  ^  ^in^in  —  ^  ^  dl0uthout 

+  N  hldiffusefr  diffuse; 


+  ^  Qht  (bulk  gas  layers) 

(19) 

Again,  the  model  simultaneously  integrates  these  equations  to  cal¬ 
culate  the  temperature  of  each  control  volume. 


2.2.3.  Heat  transfer 

Heat  is  transferred  between  layers  in  the  fuel  cell  model  by  con¬ 
duction  and  convection,  as  described  in  Section  2.1.4  by  Eqs.  (9)  and 
(10).  Heat  conduction  occurs  between  each  current  collector  and 
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the  corresponding  GDL  and  heat  convection  occurs  between  each 
bulk  gas  and  the  respective  current  collectors  and  GDLs.  Addition¬ 
ally,  convection  heat  transfer  between  the  coolant  and  the  cathode 
current  collector  removes  heat  from  the  cell. 

Convection  coefficients  are  determined  using  the  Nusselt  num¬ 
ber,  Nu : 


Nu  x  kf 


(20) 


where  kf  is  the  fluid  conduction  heat  transfer  coefficient  and  DH  is 
the  hydraulic  diameter  of  the  channel  being  modeled. 


where  Dw  is  the  diffusion  coefficient  and  t  is  the  membrane  thick¬ 
ness.  The  diffusion  is  given  as 


Dw  =  Dx  exp 


2416 


(_L 

\303 


(29) 


where 


Dx  = 


10“6 

[1  +  2(X  —  2)]  x  10-6 
[3-1 .67(7.  —  3)]  x  lO-6 
1.25  x  10“6 


for  A  <  2 
for  2  <  A  <  3 
for3  <  A  <  4.5 
for4.5  <  A 


(30) 


2.2.4.  Diffusion 

Diffusion  of  gases  from  the  bulk  gas  flow  channels  to  the  GDLs 
plays  an  important  role  in  modeling  fuel  cell  behavior  because  poor 
gas  transport  to  reaction  sites  can  lead  to  concentration  losses.  This 
model  accounts  for  diffusion  of  gases  from  the  bulk  phase  to  the 
GDL  by  calculating  a  mass  transport  coefficient: 


<§m  — 


Sh  x  Dm 
Dh 


(21) 


where  Sh  is  the  Sherwood  number  and  Dm  is  the  diffusion  coef¬ 
ficient.  This  is  then  used  to  calculate  a  total  resistance,  Rdiffuse,  to 
diffusion: 


^diffuse  —  - w  (22) 

(l/Sm)  +  (tGDL/D£f) 

where  D^f  is  the  effective  diffusion  coefficient  modified  via  the 
Bruggeman  relation: 


neff  _  pji  *  n 

um  ~  fcGDL^m 


3/2, 


(23) 


used  to  account  for  porosity  (£gdl)  and  tortuosity  in  the  GDL,  and 
tGDL  is  the  thickness  of  the  catalyst  layer.  Species  flux  is  calculated 
as 


^diffuse  —  ^diffuse  (Qz  Q)  (24) 

where  C2  and  C\  are  the  molar  concentrations  in  the  bulk  gas  phase 
and  GDL. 


2.2.5.  Water  transport 

Water  can  be  transported  through  the  ionic  membrane  via  the 
process  of  osmotic  drag  whereby  water  molecules  are  effectively 
“dragged”  from  the  anode  to  the  cathode  by  proton  transport  mech¬ 
anisms  in  the  membrane.  Additionally,  water  can  “back-diffuse” 
from  the  cathode  to  the  anode  due  to  a  concentration  gradient.  Both 
of  these  processes  are  accounted  for  in  the  current  model  because 
membrane  water  content  significantly  impacts  cell  performance 
(especially  membrane  protonic  conductivity). 

The  rate  of  water  molecules  transported  by  osmotic  drag  ($h2o) 
is  directly  proportional  to  the  number  of  protons  transferred 
through  the  membrane: 

0h2o  =  %  (25) 

where  nd  is  the  membrane  water  content  dependent  drag  coeffi¬ 
cient: 

2.5A 

"«  =  IT  (26) 

X  =  0.043  +  17.81a  -  39.85a2  +  36a3,  for0<a<l  (27) 

and  a  is  the  water  activity  as  given  by  Springer  et  al.  [19].  Water 
diffusion  due  to  the  concentration  gradient  between  the  electrodes 
is  determined  by 


2.2.6.  Electrochemistry 

The  rates  of  the  electrochemical  reactions  occurring  in  the  fuel 
cell  are  considered  instantaneous  in  the  current  model  because  the 
timescales  of  these  reactions  are  far  faster  than  those  associated 
with  the  other  physics  included  in  the  model  [  18  ].  The  electrochem¬ 
ical  potential  of  the  reactions  is  modeled  using  the  Nernst  equation 
coupled  with  equations  for  ohmic  and  activation  polarization  losses 
as  given  by  Eq.  (31),  for  anode  and  cathode  compartments  both 
operated  at  the  same  pressure: 


17  AG(T)  RT . 

V  = - ^2  +  —  In 

nF  nF 


1/2 


X, 


2  °2  pl/2 


H20 


RT 

-ac—=  In 
nF 


(c)-‘ 


-iR  (31) 


The  first  term  provides  the  theoretical  potential  based  on  Gibb’s  free 
energy  as  a  function  of  temperature.  The  change  in  Gibb’s  energy  is 
given  by 


A  G{T)= 

(-241.9495  +  0.0411  x  T+1.0641  x  10“5xT2-2.4087xl0“9xT3) 
x  1000  (32) 


The  second  term  in  Eq.  (31)  modifies  the  theoretical  potential  to 
account  for  the  operating  pressure,  temperature,  and  species  con¬ 
centrations.  The  third  term  is  a  voltage  reduction  due  to  activation 
losses  derived  from  the  Tafel  expression,  and  the  final  term  is  a 
voltage  reduction  due  to  Ohm’s  law,  where  the  resistance,  R,  is  a 
membrane  water  content  dependent  term: 


R  =  (0.005139A  -  0.00326)  exp 


1268 


(±- 

V303 


1 


273  +  T  2 


(33) 


The  rates  of  all  reactions  are  modeled  as  instantaneous  and  propor¬ 
tional  to  the  current: 


"react  =  ^  (34) 

where  i  is  the  operating  current,  Fis  Faraday’s  constant,  and  n  is  the 
number  of  moles  of  electrons  transferred  during  the  reaction  per 
mole  of  the  species  reacted.  Eq.  (34)  therefore  assumes  infinitely 
fast  reaction  with  perfect  electron  transfer  by  every  reaction. 

Concentration  polarization  has  been  neglected  in  this  model 
because  the  RFC  system  never  enters  the  operation  regime  where 
concentration  losses  become  significant.  The  maximum  current 
density  of  the  fuel  cell  is  490  mA  cm-2  when  operating  at  full  power 
of  5.5  kW. 


2.2.7.  Fuel  cell  model  parameters 

Several  parameters  used  in  this  system  model  are  given  in 
Table  3. 

2.3.  PEM  electrolyzer 


^h2o  =  A/vA 


Cc-Ca 


t 


(28) 


The  electrolyzer  model  is  based  closely  on  the  fuel  cell  model. 
It  too  is  a  bulk  model  utilizing  the  same  fundamental  equations 
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Table  3 

Fuel  cell  model  parameters 


Parameter 

Value 

Units 

Reference 

Number  of  cells  in  series 

33 

NA 

Per  design 

Exchange  current  density,  i0 

10 

Am-2 

[20] 

Activation  polarization  coefficient,  ac 

0.5 

[20] 

Operating  pressure,  P0 

1.0 

atm 

Per  design 

Dimension  of  cell 

0.2258  x  0.2258 

m 

Per  design 

Air  utilization,  /zair 

50 

% 

Per  design 

Depth  of  anode  gas  channel 

0.001 

m 

Per  design 

Depth  of  cathode  gas  channel 

0.001 

m 

Per  design 

Depth  of  cooling  channel 

0.003 

m 

Per  design 

Thickness  of  GDL,  tGDL 

0.0002 

m 

Per  design 

Thickness  of  electrolyte 

0.00015 

m 

Per  design 

Thickness  of  separator  plates 

0.002 

m 

Per  design 

GDL  porosity,  £Gdl 

0.5 

Per  design 

Separator  plate  density 

2210 

kg  m-3 

Per  design 

Separator  plate  specific  heat 

0.5 

kj(kgK)-1 

Per  design 

Electrolyte  dry  density 

2200 

kg  m-3 

Per  design 

Electrolyte  dry  equivalent  weight 

1000 

kgkmol-1 

Per  design 

Electrolyte  solid  specific  heat 

2.179 

kj(kgK)-1 

Per  design 

Separator  plate  conduction  coefficient 

0.22 

kW(mK)"1 

Per  design 

Nusselt  number  of  anode  gas,  Nu 

2.96 

[15] 

Nusselt  number  of  cathode  gas 

2.96 

[15] 

Nusselt  number  of  coolant 

7.54 

[15] 

Sherwood  number  for  anode  gas,  Sh 

12.5 

[15] 

Sherwood  number  for  cathode  gas 

12.5 

[15] 

Diffusion  coefficient  for  hydrogen 

4.1  E-5 

ms-1 

[18] 

Diffusion  coefficient  for  oxygen 

2.1E-5 

ms-1 

[18] 

Diffusion  coefficient  for  nitrogen 

2.2E-5 

ms-1 

[18] 

Diffusion  coefficient  for  hydrogen 

2.6E-5 

ms-1 

[18] 

to  calculate  temperatures,  concentrations,  and  in  this  case,  elec¬ 
trical  energy  consumption.  The  only  inputs  are  current,  ambient 
temperature,  and  water  flow  rate,  temperature,  and  pressure  to  the 
anode.  The  model  operates  on  current  control  and  calculates  elec¬ 
trolyzer  voltage  (including  losses)  and  anode  and  cathode  exit  flow 
rates,  temperatures,  and  species  compositions.  Unlike  the  fuel  cell 
model,  the  electrolyzer  does  not  include  any  separate  provisions  for 
cooling  because  the  inlet  water  is  used  to  cool  the  stack  during  all 
operating  conditions.  The  pressurized  inlet  water  is  supplied  at  a 
constant  100  g  s-1 ,  regardless  of  current,  whenever  the  electrolyzer 
is  operating.  At  maximum  power  consumption  of  3.6  kW,  the  water 
consumed  by  electrolysis  is  just  0.182  g  s_1 .  The  difference  in  water 
consumed  versus  flow  rate  provides  sufficient  cooling.  Addition¬ 
ally,  the  electrolyzer  model  has  no  anode  GDL  because  the  entire 
anode  side  of  the  cell  is  flooded.  These  differences  are  shown  in  the 
electrolyzer  model  schematic  of  Fig.  4.  It  is  assumed  that  the  anode 
bulk  flow  and  reaction  sites  are  at  the  same  temperature,  and  that 
as  long  as  the  water  flow  rate  to  the  electrolyzer  is  sufficient,  then 


Anode  Current  Collector 
Anode  Bulk  Water/O; 


Membrane 
Cathode  GDL 


Cathode  Bulk  Gas 
Cathode  Current  Collector 


Fig.  4.  Diagram  showing  six  control  volumes  modeled  for  the  electrolyzer. 


Table  4 

Electrolyzer  model  parameters 

Parameter  Value  Units  Reference 


Number  of  cells  in  series 
Exchange  current  density,  i0 
Activation  polarization  coefficient,  ac 
Operating  pressure,  P0 
Dimension  of  cell 
Air  utilization,  /zair 
Depth  of  anode  gas  channel 
Depth  of  cathode  gas  channel 
Thickness  of  GDL,  tGDL 
Thickness  of  electrolyte 
Thickness  of  separator  plates 
GDL  porosity,  £Gdl 
Separator  plate  density 
Separator  plate-specific  heat 
Electrolyte  dry  density 
Electrolyte  dry  equivalent  weight 
Electrolyte  solid-specific  heat 
Separator  plate  conduction  coefficient 
Nusselt  number  of  anode  fluid,  Nu 
Nusselt  number  of  cathode  gas 
Sherwood  number  for  cathode  gas 
Diffusion  coefficient  for  hydrogen 
Diffusion  coefficient  for  oxygen 
Diffusion  coefficient  for  hydrogen 


15 

NA 

Per  design 

0.01 

Am-2 

[20] 

0.5 

[20] 

6.8 

atm 

Per  design 

0.13  x  0.13 

m 

Per  design 

50 

% 

Per  design 

0.01 

m 

Per  design 

0.001 

m 

Per  design 

0.0002 

m 

Per  design 

0.00015 

m 

Per  design 

0.002 

m 

Per  design 

0.5 

Per  design 

2210 

kg  nrr3 

Per  design 

0.5 

kJ(kgK)-1 

Per  design 

2200 

kg  ITT3 

Per  design 

1000 

kg  kmol 

Per  design 

2.179 

kJ(kgK)-1 

Per  design 

0.22 

kW(mK)"1 

Per  design 

7.54 

[15] 

2.96 

[15] 

12.5 

[15] 

4.1  E-5 

ms-1 

[18] 

2.1E-5 

ms-1 

[18] 

2.6E-5 

ms-1 

[18] 

there  is  always  liquid  water  at  the  anode  catalyst  sites  with  no  diffu¬ 
sion  limitations  in  this  region.  The  MEA  therefore  consists  of  only  a 
membrane  and  the  cathode  GDL.  An  additional  difference  between 
the  fuel  cell  and  electrolyzer  models  is  that  no  nitrogen  is  included 
in  any  of  the  electrolyzer  control  volumes.  Otherwise,  the  same  set 
of  equations  as  those  used  in  the  fuel  cell  model  are  simultaneously 
solved  in  each  time  step  in  the  electrolyzer  model. 


2.3.1.  Electrolyzer  electrochemistry 

As  in  the  fuel  cell  component  model,  the  rates  of  the  electro¬ 
chemical  reactions  in  the  electrolyzer  are  assumed  instantaneous 
because  the  timescales  of  these  reactions  are  far  faster  than  those 
associated  with  the  other  physics  in  the  model.  The  electrochemical 
potential  of  the  electrolyzer  is  modeled  using  the  Nernst  equa¬ 
tion  coupled  with  equations  for  ohmic  and  activation  polarization 
losses: 


17  AG(T)  RT . 

V  = - ^  +  -=  In 

nF  nF 


XH,Xo/2  , 

M2  02  pi/2 


X, 


h20 


RT.  ,  . 
+  ac-ln(^ 


tt)+“ 


(35) 


The  only  differences  between  Eq.  (31)  for  the  fuel  cell  and  Eq.  (35) 
for  the  electrolyzer  are  sign  changes  on  the  activation  and  ohmic 
loss  terms.  The  rates  of  all  reactions  are  again  proportional  to  the 
current: 

bread:  =  (36) 

Concentration  polarization  has  been  neglected  due  to  the  fact 
that  the  electrolyzer  never  operates  at  high  current  density  con¬ 
ditions.  Ohmic  resistance  in  the  electrolyzer  is  a  function  of  water 
content.  The  electrolyzer  model  ends  up  predicting  ohmic  losses 
that  are  relatively  low  because  the  anode  is  always  flooded,  reduc¬ 
ing  water  content  variation  in  the  membrane.  Like  the  fuel  cell, 
both  anode  and  cathode  compartments  are  pressurized  at  the  same 
level  during  operation.  By  pressurizing  the  anode  inlet  water,  the 
use  of  a  separate  hydrogen  compressor  that  would  otherwise  be 
needed  to  generate  the  pressure  required  by  the  hydride  tank  can 
be  sidestepped. 


2.3.2.  Model  parameters 

Table  4  lists  various  parameters  used  in  the  electrolyzer  model. 
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Coolant  Out 


Line  of  axial  symmetry 


Tube  wall 


Coolant  In 


Fig.  5.  Heat  transfer  pathways  of  radiator  module. 


2.4.  Radiator 

A  first  principles  model  of  a  simple,  one  dimensional,  tube  and 
fin  radiator  consisting  of  four  adjacent  nodes  has  been  developed 
to  simulate  the  system  radiator.  This  radiator  serves  to  cool  the 
hydride  tank  during  charging  and  to  cool  the  fuel  cell  during  power 
mode.  Fig.  5  shows  a  diagram  of  the  heat  transfer  modes  in  the 
tube  of  the  radiator.  Convection  heat  transfer  occurs  between  the 
coolant  and  the  inner  tube  wall,  heat  conduction  occurs  between 
the  tube  nodes,  and  convective  heat  transfer  with  the  environment 
occurs  on  the  outer  tube  wall.  The  gas  exit  temperatures  of  each 
node  are  used  to  calculate  the  temperature  differences  needed  for 
the  convection  heat  transfer  equation  (37)  between  the  coolant  and 
the  tube  wall.  The  equation  for  heat  convection,  Qconvection*  between 
coolant  and  tube  at  each  node  is  calculated  as 

Qconvection  =  ^coolantA(^coolant_out  -Ttube)  (37) 

where  hcoolant  is  the  coolant  convection  heat  transfer  coefficient, 
A[  is  the  internal  tube  surface  area,  and  Tcoolant  out  and  Ttube  are  the 
exit  coolant  temperature  and  tube  nodal  temperature,  respectively. 
Convection  heat  transfer  between  the  tube  and  the  environment  is 

Qconvection  =  ^environ71e(Ftube  —  ^environ)  (38) 

where  henv iron  is  the  ambient  environmental  convection  heat  trans¬ 
fer  coefficient,  Ae  is  the  external  tube  surface  area,  and  Ttube  and 
^environ  are  the  tube  nodal  temperature  and  ambient  temperature, 
respectively.  Conduction  along  the  length  of  the  tube  is  calculated 
similarly  to  Eq.  (9).  Table  5  lists  the  parameter  values  used  in  the 
radiator  model. 

2.5.  Ancillary  components 

Several  other  important  components  have  been  modeled  empir¬ 
ically  to  complete  the  system.  A  simple  on-off  100  W  radiator  fan 
provides  the  forced  convection  necessary  to  remove  heat  from  the 
radiator.  When  the  fan  is  off,  the  radiator  is  cooled  by  natural  con- 


Table  5 

Radiator  model  parameters 


Parameter 

Value 

Units 

Tube  conduction 
coefficient,  k 

237 

W(mK)-1 

Coolant  to  tube  convection 
coefficient,  /ict 

500 

W(m2  K)-1 

Air  to  tube  convection 
coefficient,  hat 

50  if  fan  is  on;  0.5  if  fan  is  off 

W(m2  K)_1 

Tube  length,  lt 

2.54 

m 

Tube  thickness,  tt 

0.001 

m 

Tube  inner  diameter 

0.024 

m 

Tube  density 

2700 

kg  ITT  3 

Tube  specific  heat 

903 

J(kgK)-1 

Coolant  density 

1000 

kg  ITT3 

Coolant-specific  heat 

4200 

J(kgK)-1 

Coolant  flow  rate 

0.33 

kgs-1 

vection  with  the  ambient  environment.  When  the  fan  is  turned  on, 
the  convection  heat  transfer  coefficient  between  the  radiator  tube 
and  the  environment  is  increased  to  simulate  forced  convection. 
Another  ancillary  component  is  a  constant  power  200  W  coolant 
pump  that  circulates  the  coolant  fluid  at  a  continuous  0.33 kgs-1 
during  all  modes  of  system  operation.  The  final  ancillary  compo¬ 
nent  modeled  is  a  variable  speed  fuel  cell  inlet  air  blower.  The 
blower  power  is  proportional  to  the  cube  of  the  required  air  flow 
[21]  with  a  maximum  power  level  of  540  W  at  an  air  flow  rate  of 
0.218  mol  s-1.  The  required  air  flow  rate  at  any  fuel  cell  operating 
point  is  determined  from  the  fuel  cell  current  level  and  the  desired 
oxidant  utilization  factor,  /xair: 


hair  — 


i 

4fytairX02 


(39) 


Xo2  is  the  oxygen  mole  fraction  in  air  and  the  4  is  due  to  4  mol  of 
electrons  transferred  per  mole  of  oxygen  gas. 


3.  System  interfaces 

Fig.  6  shows  a  diagram  of  the  physical  interfaces  between  RFC 
components  in  the  current  model.  During  charging  mode,  the  elec¬ 
trolyzer  is  supplied  with  electrical  current  and  deionized  water 
while  the  hydrogen  produced  is  absorbed  exothermically  by  the 
hydride  tank.  The  heat  from  this  reaction  is  removed  by  a  circulating 
water  bath  cooled  by  the  radiator.  In  power  mode,  the  electrolyzer 
is  turned  off  and  the  fuel  cell  produces  power  from  hydrogen  sup¬ 
plied  by  the  endothermic  desorption  of  hydrogen  from  the  metal 
hydride  bed.  Simultaneously,  the  heat  generated  by  the  fuel  cell  is 
removed  via  liquid  cooling  and  fed  to  the  hydride  tank  water  bath  to 
simultaneously  provide  the  necessary  energy  for  hydrogen  evolu¬ 
tion  and  be  cooled  by  this  endothermic  process.  The  exiting  water 
is  then  further  cooled  by  the  radiator  before  re-entering  the  fuel 
cell. 

The  three  numbered  boxes  in  Fig.  6  represent  model  processes  in 
which  the  dynamics  of  the  extant  physics  are  not  accounted  for.  Box 
number  1  represents  an  enthalpy  exchange  device  (e.g.,  enthalpy 
wheel)  that  is  simplified  to  set  the  temperature  and  humidity  of 
the  incoming  air  to  65  °C  and  30%  relative  humidity,  respectively, 
regardless  of  the  ambient  conditions.  Box  number  2  simulates  a 
throttling  and  enthalpy  exchange  process  wherein  hydrogen  pres¬ 
sure  drops  from  the  internal  hydride  tank  pressure  down  to  the  fuel 
cell  inlet  pressure  (110  kPa),  and  the  temperature  and  humidity  are 
increased  to  65  °C  and  75%  relative  humidity,  respectively,  via  an 
assumed  enthalpy  exchange  with  the  fuel  cell  cathode  exhaust. 
The  total  enthalpy  of  the  stream  exiting  box  number  2  is  equal 
to  the  total  enthalpy  of  the  stream  entering  box  number  2.  The 
internal  hydride  tank  pressure  is  governed  by  the  tank  tempera¬ 
ture  (Eq.  (2))  and  is  always  controlled  to  be  greater  than,  or  equal 
to  the  fuel  cell  inlet  pressure,  as  explained  in  the  following  controls 
section.  The  maximum  tank  pressure  is  controlled  to  be  689.5  kPa. 
Therefore,  the  largest  possible  pressure  drop  during  the  throttling 
expansion  would  be  579.5  kPa.  Using  a  Joule-Thomson  coefficient 
for  room  temperature  hydrogen  of  4.0E-4KkPa-1  [22],  the  tem¬ 
perature  change  during  the  expansion  due  to  the  Joule-Thomson 
effect  can  be  shown  to  be  negligible,  which  supports  the  current 
simplification  of  this  process  as  follows  [23]: 

„.(S), 

AT  =  /MjAP  =  4  X  10"4  KkPa-1  X  579.5  kPa  =  0.23 1<  (41 ) 

Box  number  3  simulates  a  similar  reduction  in  pressure  and  a 
reduction  of  humidity  process  wherein  the  outlet  electrolyzer  pres¬ 
sure  (689  kPa)  is  reduced  to  the  necessary  tank  inlet  pressure,  and 
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Fig.  6.  Diagram  of  RFC  system  model.  Numbered  boxes  represent  state  changes  in  which  the  physics  are  not  fully  modeled. 


the  gas  is  dried  to  0%  relative  humidity  (e.g.,  by  using  a  desiccant). 
Again,  the  electrolyzer  pressure  is  controlled  to  always  be  greater 
than,  or  equal  to  the  tank  inlet  pressure. 

4.  Results 


4.1.  System  control 


Several  interesting  situations  can  arise  during  operation  of  this 
system.  During  the  initial  stages  of  system  charging,  the  elec¬ 
trolyzer  maximum  hydrogen  production  rate  will  be  the  bottleneck 
because  the  hydride  tank  is  cool,  empty,  and  can  accept  a  larger 
flow  of  hydrogen  than  the  electrolyzer  can  produce.  However,  as 
the  hydride  tank  fills  and  heats  during  charging,  the  available  fill 
rate  slows,  eventually  reaching  the  point  where  the  rate  at  which 
the  tank  can  fill  drops  below  the  rate  at  which  the  electrolyzer  can 
produce  hydrogen.  The  model  accounts  for  this  by  calculating  a  line 
pressure  between  the  electrolyzer  and  the  hydride  tank  using  the 
Ideal  Gas  Law: 


P  = 


nRT 


(42) 


where  the  number  of  moles,  n,  in  the  line  volume  can  be  found  by 
integrating  a  molar  balance: 


dn 

=  ne\ec  —  ntank 


(43) 


Here  nelec  is  the  molar  flow  rate  of  hydrogen  produced  by  the 
electrolyzer  and  ntank  is  the  molar  flow  rate  of  hydrogen  enter¬ 
ing  the  hydride  tank.  The  electrolyzer  always  produces  hydrogen 
output  at  689  kPa  because  the  water  supplied  to  the  electrolyzer 
is  pressurized  to  689  kPa.  Initially,  the  line  pressure  stays  well 
below  689  kPa  because  the  rate  of  hydrogen  absorbed  by  the  tank 


is  nearly  identical  to  the  rate  produced  by  the  electrolyzer,  and  lit¬ 
tle  mass  can  accumulate  in  the  line  between  the  two  components. 
During  this  operating  condition,  the  tank  fills  at  a  linear  rate  pro¬ 
portional  to  the  electrolyzer’s  maximum  output.  As  the  tank  fill  rate 
slows  due  to  tank  heating  and  fill  level  (Eqs.  (1)  and  (2)),  the  line 
pressure  builds,  eventually  reaching  689  kPa.  As  this  threshold  is 
approached,  a  proportional  and  triggered  integral  controller  based 
on  pressure  feedback  reduces  the  current  powering  the  electrolyzer 
in  order  to  maintain  the  maximum  line  pressure.  The  system  tran¬ 
sient  response  to  this  charging  situation  is  shown  in  Fig.  7.  Note 
that  electrolyzer  current  is  reduced  considerably  once  the  pressure 
reaches  7  bar,  but,  that  the  hydride  tank  continues  to  fill  (albeit  at 
a  slower  rate). 
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Fig.  7.  RFC  transient  response  during  charging  operation  with  increasing  tank  fill 
level. 
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Fig.  8.  RFC  temperature  responses  during  charging  operation  with  increasing  tank 
fill  level. 


Fig.  8  shows  temperature  information  during  the  same  filling 
event.  The  temperature  profiles  for  the  electrolyzer  membrane  and 
inner  node  of  the  hydride  tank  are  plotted  on  the  left  axis  and 
the  temperature  profiles  for  the  surface  of  the  hydride  tank  and 
the  coolant  water  into  and  out  of  the  hydride  cooling  chamber 
are  plotted  on  the  right,  to  better  distinguish  the  relatively  small 
temperature  differences  between  them.  It  can  be  seen  that  the  elec¬ 
trolyzer  reaches  an  operating  temperature  of  roughly  80  °C  within 
several  minutes,  and  maintains  this  temperature  throughout  the 
fill  event.  The  inner  node  of  the  hydride  tank  warms  fairly  steadily 
until  this  portion  of  the  tank  becomes  saturated  (around  100  min) 
at  which  time  the  temperature  steadily  decreases  because  the 
exothermic  reaction  has  stopped.  Despite  the  inner  tank  tempera¬ 
ture  rising  by  almost  40  °C,  the  tank  surface  temperature  only  rises 
6.5  °C  due  to  contact  with  the  coolant  bath  and  the  poor  heat  trans¬ 
fer  properties  of  the  powdered  metal  hydride  material.  The  coolant 
temperatures  show  that  relatively  little  energy  is  transferred  from 
the  tank  to  the  coolant  (because  of  the  small  temperature  difference 
between  the  coolant  inlet  and  outlet  temperature),  underscoring 
the  rate-determining  effect  of  hydride  tank  heat  transfer. 

Fig.  9  shows  a  diagram  of  the  strategy  for  controlling  electrolyzer 
current.  This  controller  is  designed  to  complete  the  charging  pro¬ 
cess  as  quickly  as  possible  by  first  operating  at  the  maximum 
electrolyzer  output,  and  then  transitioning  to  operate  at  the  maxi¬ 
mum  allowable  system/line  pressure  as  the  hydride  tank  filling  rate 
becomes  the  system  bottleneck.  This  charging  strategy  is  logical  for 
an  onboard  military  applications  where  alternator  power  is  plenti¬ 
ful  and  maximum  charge  rate  is  more  important  than  overall  system 
efficiency.  However,  one  can  easily  imagine  different  charging  pri- 


maximum  allowable  pressure 

Fig.  9.  Diagram  of  electrolyzer  current  control  strategy. 
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Fig.  10.  Diagram  of  fuel  cell  current  and  gas  flow  rate  control  strategy. 

orities  (e.g.,  to  maximize  efficiency  or  minimize  cost  by  reducing 
component  sizes)  not  considered  here. 

Another  intriguing  situation  can  occur  during  fuel  cell  operation 
if  the  power  requested  from  the  fuel  cell  is  greater  than  that  avail¬ 
able  from  the  tank.  An  approach  similar  to  that  of  the  electrolyzer  is 
taken  whereby  the  gas  pressure  in  a  line  volume  between  the  tank 
outlet  and  fuel  cell  inlet  is  simulated  based  on  the  Ideal  Gas  Law 
and  a  dynamic  molar  balance  equation.  The  fuel  cell  inlet  pressure 
is  110  kPa  requiring  that  the  line  pressure  not  drop  below  this  value. 
Additionally,  fuel  utilization  is  controlled  to  maintain  cell  voltage 
and  ensure  efficient  operation. 

The  fuel  cell  control  strategy  uses  pressure  feedback  to  control 
current.  If  the  line  pressure  is  greater  than  the  minimum  pressure, 
the  current  will  default  to  the  desired  current.  If  the  tank  cannot 
supply  enough  fuel,  the  line  pressure  will  drop,  and  the  current  will 
be  reduced  proportionally  to  match  the  hydrogen  supply  available. 

The  secondary  utilization  control  loop  uses  a  utilization  error 
signal  to  proportionally  control  the  fuel  flow  rate.  This  double 
controller  configuration  enables  current  and  utilization  to  be  con¬ 
trolled  even  as  the  hydrogen  supply  is  diminished.  Fig.  10  shows  a 
diagram  illustrating  both  fuel  cell  controllers. 

Fig.  11  shows  the  system  response  during  low  fuel  level  opera¬ 
tion.  The  current  is  stable  at  the  desired  level  of  100  A  until  the  line 
pressure  drops  to  110  kPa,  at  which  point  the  current  and  hydrogen 
flow  rate  drop  to  ensure  sufficient  pressure  and  to  maintain  80%  fuel 
utilization.  The  plot  demonstrates  how  excess  hydrogen  pressure  is 
initially  available  as  the  fuel  cell  waste  heat  aides  the  hydride  des¬ 
orption  reaction.  The  hydride  tank  empties  at  a  nearly  linear  rate  as 
the  system  operates  at  constant  power.  However,  when  the  tank  fill 
level  falls,  the  desorption  reaction  slows  and  the  pressure  available 
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Fig.  11.  RFC  transient  response  during  power  operation  with  diminishing  fuel  level 
in  the  hydride  tank. 
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level. 


Fig.  13.  System  round-trip  efficiency  as  it  varies  with  system  power  output. 


from  the  tank  drops.  Note  that  a  nearly  constant  fuel  utilization 
of  80%  is  maintained  throughout  the  operation.  The  inconsistent 
line  pressure  decrease  during  the  first  20  min  of  operation  can  be 
accounted  for  by  a  drastic  initial  pressure  decrease  as  the  fuel  cell 
begins  using  fuel,  then  a  leveling  off  portion  as  the  waste  heat  from 
the  fuel  cell  begins  to  aid  hydrogen  flow  from  the  tank,  then  a  sec¬ 
ond,  less  severe  pressure  reduction  caused  by  decreased  hydrogen 
supply  as  the  tank  level  drops. 

Fig.  12  shows  temperature  results  for  the  fuel  cell  power  event. 
The  temperatures  of  the  fuel  cell  membrane,  hydride  tank  surface, 
coolant  into  the  tank,  and  coolant  out  of  the  tank  all  follow  the  same 
profile.  They  rise  steadily  for  the  first  20  min,  then  begin  to  level  off 
and  even  slightly  decrease  after  the  fuel  cell  current  is  reduced  by 
the  controller.  The  internal  temperature  of  the  hydride  tank  falls 
below  room  temperature  as  the  endothermic  hydrogen  desorption 
reaction  occurs,  but  cannot  fully  utilize  the  fuel  cell  waste  heat  due, 
again,  to  the  poor  heat  transfer  characteristics  of  the  tank.  The  radi¬ 
ator  plays  no  roll  in  this  scenario  because  the  fuel  cell  never  reaches 
its  design  set  point  temperature  of  80  °C. 

4.2.  System  efficiencies 

The  RFC  system  presented  herein  is  relatively  simple  and  self- 
contained.  Apart  from  the  electrolyzer,  the  only  energy  consuming 
devices  are  a  variable  power  0-540  W  blower  for  fuel  cell  inlet  air, 
a  100  W  radiator  fan,  and  a  200  W  water  pump.  System  efficiency 
for  any  operating  period  can  therefore  be  defined  by  Eq.  (44): 

_  £fuel  cell  —  ^blower  —  ^FC_pump  —  ^FCJan  r  a  a\ 

^electrical  r  ,  r  ,  n  (44) 

^electrolyzer  ^E.pump  ~r  cE_fan 

where  Efueiceii  and  ^electrolyzer are  the  quantities  of  electrical  energy 
supplied  by  the  fuel  cell  and  delivered  to  the  electrolyzer,  respec¬ 
tively.  The  electrical  energy  required  for  the  blower,  Flower*  is 
accounted  for  as  well  as  the  energy  used  by  the  fan  and  pump  dur¬ 
ing  electrolysis,  ^EJan*  ^E_pump»  and  during  fuel  cell  operation,  EFC  fan, 

^FC.pump- 

System  round-trip  efficiencies  have  been  tabulated  for  various 
system  power  levels  as  shown  in  Fig.  13.  Due  to  the  diminishing 
hydrogen  desorption  rate  as  the  level  in  the  hydride  tank  falls,  no 
power  set  point  can  be  maintained  for  an  entire  tank  of  fuel.  For 
example,  a  system  power  output  of  3.0  kW  can  be  maintained  as 
the  tank  drains  from  100%  fill  level  to  50%,  or  a  system  power  of  just 
50  W  can  be  maintained  all  the  way  down  to  a  tank  level  of  5.5%,  but 
once  each  of  these  respective  levels  is  reached,  the  hydrogen  flow 
rate  becomes  so  small  that  the  fuel  cell  current  controller  must 


reduce  the  current  output  to  maintain  the  necessary  fuel  line  pres¬ 
sure.  To  compare  system  efficiencies  at  different  power  levels  on  the 
same  basis,  the  efficiencies  shown  in  Fig.  13  are  calculated  by  oper¬ 
ating  at  constant  power  from  a  100%  full  tank  only  until  that  partic¬ 
ular  power  level  can  no  longer  be  maintained,  and  then  switching 
to  charge  mode  to  fill  the  tank  back  to  100%.  This  means  that  the 
total  energy  (power  multiplied  by  time)  stored  and  delivered  for 
each  point  on  the  efficiency  curve  is  not  necessarily  the  same. 

The  lower  efficiencies  observed  at  low  power  levels  are  due  to 
relatively  high  parasitic  loads,  such  as  the  constant  200  W  cooling 
pump  load.  The  decrease  in  efficiency  at  high  power  levels  is  due 
to  reduced  fuel  cell  and  electrolyzer  efficiency  at  higher  loads  com¬ 
bined  with  increased  parasitic  losses  from  the  variable  power  fuel 
cell  air  inlet  blower.  It  is  important  to  note  that  the  cooling  sys¬ 
tem  fan  never  operated  during  any  of  the  power  cycles  presented 
in  Fig.  13  because  the  fuel  cell  never  reached  the  design  operat¬ 
ing  temperature  of  80  °C.  Though  the  internal  fuel  cell  temperature 
gradually  increased  throughout  operation  at  each  power  level  (in 
a  similar  fashion  to  that  shown  in  Fig.  12),  the  endothermic  hydro¬ 
gen  desorption  reaction  in  the  hydride  tank  removed  sufficient  heat 
from  the  coolant  stream  to  keep  the  fuel  cell  temperature  below  its 
high  temperature  limit.  A  larger  hydrogen  storage  tank  may  have 
allowed  the  simulated  operation  to  continue  long  enough  for  the 
fuel  cell  to  eventually  reach  the  maximum  operating  temperature. 
This  cooling  phenomenon  can  help  explain  the  fairly  flat  mid-range 
efficiency  shown  in  Fig.  13.  In  the  current  PEM  fuel  cell  model,  effi¬ 
ciency  decreases  as  power  level  increases  and  fuel  cell  efficiency 
increases  with  increasing  temperature.  Therefore,  since  the  fuel  cell 
temperature  never  reaches  an  upper  boundary  during  the  simula¬ 
tions,  the  efficiency  losses  due  to  higher  current  density  operation 
are  offset  by  the  efficiency  gains  reached  through  the  gradually 
higher  temperature  operation. 

Figs.  14  and  15  show  the  power  outputs  and  inputs  for 
each  of  the  system  components  during  transient  operation  of  a 
charge/discharge  cycle  from  an  initial  hydride  tank  level  of  25-95% 
and  back  to  25%  at  a  constant  current  set  point.  Both  plots  are  iden¬ 
tical  during  the  tank  filling  portion  of  the  operation  because  the 
electrolyzer  fills  the  tank  from  25  to  95%  at  the  maximum  rate  as 
described  previously.  In  Fig.  14  (Case  1 ),  after  reaching  a  tank  level 
of  95%,  the  electrolyzer  is  turned  off  and  the  fuel  cell  is  operated  at 
approximately  1.35  kW  (50  A)  until  the  tank  has  been  drained  back 
to  a  25%  level.  After  subtracting  parasitic  pump  and  fan  loads,  this 
equates  to  a  system  output  power  of  roughly  1.15  kW.  Fig.  14  shows 
that  the  system  is  capable  of  maintaining  this  power  output  level 
for  the  duration  of  the  test  (275  min). 
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Fig.  14.  Transient  RFC  system  response  (Case  1)  showing  component  power  levels 
beginning  with  charging  mode  and  switching  to  power  mode  after  377  min.  The  sys¬ 
tem  then  maintains  a  constant  power  output  of  roughly  1.15  kW  until  the  completion 
of  the  simulation. 


Fig.  16.  Maximum  theoretical  volumetric  and  gravimetric  power  densities  for  an 
RFC  system  utilizing  components  that  possess  performance  characteristics  similar 
to  those  commercially  available. 
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Fig.  15.  Transient  RFC  system  response  (Case  2)  showing  component  power  levels 
beginning  with  charging  mode  and  switching  to  power  mode  after  377  min.  The 
system  then  attempts  to  maintain  a  constant  power  output  of  2.4  kW. 


Fig.  15  shows  Case  2  which  has  a  similar  charge  and  discharge 
cycle  but  with  a  larger  fuel  cell  power  output  of  approximately 
2.70  kW  ( 100  A)  corresponding  to  a  system  power  of  approximately 
2.4  kW.  It  can  be  seen  that  the  system  is  unable  to  maintain  this 
higher  power  level  for  the  duration  of  the  test.  Fig.  15  shows  a 
moderate  output  power  improvement  as  the  fuel  cell  performance 
increases  due  to  the  steady  temperature  rise  of  the  stack,  but  after 
107  min,  the  current  controller  begins  to  reduce  fuel  cell  power  as 
hydrogen  flow  becomes  limited.  At  this  point,  the  fuel  cell  waste 
heat  carried  by  the  cooling  fluid  to  the  hydride  tank  is  insufficient 
to  remove  the  quantity  of  gas  required  by  the  fuel  cell.  This  occurs 
because  the  hydrogen  flow  rate  from  the  tank  slows  as  the  tank  level 
drops.  It  should  also  be  noted  that  the  radiator  cooling  fan  was  not 
powered  because  the  coolant  temperature  never  reached  the  “fan 
on”  threshold. 

Table  6  lists  this  energy  data  and  the  corresponding  overall  sys¬ 
tem  efficiencies  for  both  simulation  cases. 


The  efficiency  of  29.8%  calculated  in  Table  6  for  operation  at 
1.15  kW  is  slightly  higher  than  the  efficiency  predicted  by  Fig.  13 
for  the  same  power  output  (29.1%).  This  is  because  the  charge  por¬ 
tion  of  the  cycle  evaluated  in  Fig.  13  fills  the  tank  all  the  way  to  100%, 
instead  of  just  95%.  This  small  difference  in  fill  level  has  a  significant 
effect  on  overall  efficiency  because  of  the  constant  parasitic  load  of 
the  200  W  water  pump.  Although  the  electrolyzer  is  operating  very 
efficiently  during  the  final  5%  of  the  fill  because  the  hydride  tank 
absorption  rate  at  this  point  is  so  slow  (low  current  density  opera¬ 
tion  of  the  electrolyzer),  the  water  pump  continues  to  work  at  full 
power,  lowering  the  system  efficiency.  A  similar  effect  occurs  with 
Case  2  for  operation  at  2.40  kW.  These  types  of  insights,  often  not 
included  in  systems  analyses  that  only  simulate  the  major  compo¬ 
nents,  are  important  to  understand  for  the  design  and  operation  of 
complete  RFC  systems. 

4.3.  Energy  and  power  density 

The  RFC  system  studied  in  this  analysis  could  be  built  by 
using  currently  available  commercial  products.  A  4  kW  Hydrogen- 
ics  HyLYZER  PEM  electrolyzer  and  a  5  kW  Flydrogenics  FlyPM  fuel 
cell  [24]  could  be  used  in  conjunction  with  10,000  standard  liters  of 
hydrogen  storage  in  the  form  of  three  metal  hydride  tanks  available 
from  GfE  Metalle  und  Materialien,  GmbH  of  Ntirnberg,  Germany. 
The  dimensions  of  each  component,  as  well  as  those  for  the  stan¬ 
dard  military  lead-acid  battery  (“6TMF”)  currently  used  in  the 
Stryker  vehicle  [25],  and  a  commercially  available  lead-acid  battery 
(“Hawker  Genesis  G70EP”)  [26]  are  given  in  Table  7.  The  Hawker 
Genesis  is  included  because  of  the  wealth  of  performance  data 
available  from  the  manufacturer  and  in  the  literature.  Numbers  in 
the  gray  regions  are  actual  values  provided  by  manufacturers.  Val¬ 
ues  listed  in  the  non-shaded  areas  are  estimates  used  in  the  current 
model. 

The  gravimetric  and  volumetric  energy  densities  given  in  Table  7 
clearly  indicate  that  the  hypothetical  RFC  system  analyzed  herein 
is  not  currently  competitive  with  lead-acid  battery  technology  on  a 
pure  energy  density  basis.  Likewise,  an  RFC  system  cannot  produce 
nearly  as  much  power  on  a  mass  or  volume  basis  as  a  lead-acid  bat- 


Table  6 

Energies  and  overall  system  efficiencies  for  the  two  RFC  operation  cases 


Operating  condition 

^electrolyzer  (MJ) 

^E.pump(MJ) 

£E.fan(MJ) 

^fuelcell  (MJ) 

^blower  (MJ) 

^FC.pump  (MJ) 

^FCfan  (MJ) 

^electrical  (^>) 

Case  1 

57.25 

4.536 

2.268 

22.430 

0.071 

3.276 

0.0 

29.8 

Case  2 

57.25 

4.536 

2.268 

21.893 

0.2507 

1.848 

0.0 

30.9 
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Table  7 

RFC  system  component  and  lead-acid  battery  dimensions,  and  calculated  energy  and  power  densities 


Component 

Volume 

Mass 

Energy 

Power 

Volumetric 

Gravimetric 

Power 

Power 

(Liter) 

(kg) 

Output 

Output 

Energy 

Energy 

/Vol 

/Mass 

(MJ) 

(kW) 

Density 

Density 

(kW  L-1) 

(kW  kg-i) 

(kJ  L'1) 

(kJ  kg-0 

Lead  acid  battery,  6TMF 

15 

33.5 

5.18a 

NA 

346 

156 

NA 

NA 

Lead  acid  battery,  G70EP 

9.8 

24.3 

3.38a 

3.68b 

345 

139 

0.38 

0.14 

Hydrogenics  HyPM10c 

78 

44.0 

Hydrogenics  HyLYSER  2.0 

34 

30.0 

10,000  liter  hydride  tank 

20 

125 

Tank  cooling  bath  and  enclosure 

25 

10.0 

System  total 

137 

209 

18. 9d 

5.0 

138 

90 

0.036 

0.024 

aEnergy  output  during  C20  discharge  cycle.  bPower  output  sustainable  for  5  min.  cProduct  specifications  were  available  only  for  a  10  kW  fuel  cell,  not  a  5  kW  cell  as  studied 
in  this  analysis.  Therefore,  the  size  and  weight  of  the  10  kW  fuel  cell  were  halved  to  estimate  the  specifications  for  the  smaller  cell.  dSystem  energy  output  during  continuous 
1.15  kW  operation. 


Fig.  17.  Maximum  theoretical  volumetric  and  gravimetric  energy  densities  for  an 
RFC  system  utilizing  components  that  possess  performance  characteristics  similar 
to  those  commercially  available. 

tery.  However,  the  storage  capability  of  the  RFC  is  close  enough  to 
that  of  the  battery  that  one  can  envision  niche  applications  where 
the  theoretically  superior  cyclic  durability,  modular  design,  and  lack 
of  self-discharge  properties  of  the  RFC  system  may  outweigh  the 
energy  density  deficiency  as  compared  to  traditional  battery  tech¬ 
nology.  The  scope  of  such  niche  applications  will  only  increase  as 
RFC  components  and  systems  are  optimized.  In  addition,  improve¬ 
ments  in  the  energy  and  power  density  of  RFC  system  components 
are  expected  in  the  future,  which  could  significantly  improve  com¬ 
petitiveness  with  battery  systems. 

Because  of  the  disassociated  nature  of  the  power  and  energy  per¬ 
formance  characteristics  of  an  RFC  system,  it  is  easy  to  project  the 
theoretical  maximum  energy  and  power  densities  attainable  with 
given  component  technology  performance.  The  energy  density 
can  clearly  never  exceed  the  basic  energy  density  of  the  particu¬ 
lar  hydrogen  storage  tank  used.  Likewise,  the  power  output  can 
never  exceed  the  power  output  of  the  fuel  cell.  A  limiting  case  for 
both  scenarios  is  reached  as  the  tank  or  fuel  cell  size  is  increased 
to  the  point  where  the  sizes  of  the  other  components  become 
insignificant,  as  shown  in  Figs.  16  and  17.  These  projections  assume 
that  the  power  output  and  hydrogen  stored  increase  linearly  with 
size  and  mass  of  the  fuel  cell  and  tank,  respectively.  Though  fuel 
cell  power  does  increase  linearly  with  electrode  surface  area,  and 
hydrogen  storage  does  increase  linearly  with  the  mass  of  metal 
hydride,  an  actual  larger  system  may  offer  greater  performance 
gains  than  predicted  here  due  to  the  non-linear  performance-to- 


size  changes  of  balance  of  plant  equipment  components,  such  as 
blowers  and  pumps.  These  projections  can  therefore  be  consid¬ 
ered  a  “worst  case”  scenario.  Figs.  16  and  17  show  maximum  power 
densities  approaching  0.1  kWkg-1  and  0.06  kWL-1  and  maximum 
energy  densities  of  more  than  800  kJL-1  and  150kjkg_1,  respec¬ 
tively.  This  shows  that  the  potential  volumetric  energy  density  of  an 
RFC  system  far  exceeds  that  of  a  lead-acid  battery  while  the  poten¬ 
tial  gravimetric  energy  density  is  roughly  equivalent,  with  current 
metal  hydride,  PEM  fuel  cell  and  PEM  electrolyzer  performance 
characteristics. 

5.  Summary  and  conclusions 

An  RFC  system  model  has  been  built  using  first  principles 
approaches  for  all  of  the  major  system  components  and  simple 
feedback  control  strategies  for  the  system  interfaces.  The  model 
has  been  exercised  and  shown  to  demonstrate  reasonable  system 
transient  responses  during  all  modes  of  operation.  The  following 
conclusions  can  be  drawn  from  simulations  conducted  with  this 
model: 

•  Charge  rate  is  at  first  governed  by  electrolyzer  power.  The  hydride 
tank  fill  rate  will  then  slow  to  a  point  where  it  becomes  the  rate 
limiting  factor.  The  fill  rate  of  the  first  portion  of  the  charging  bot¬ 
tleneck  can  be  increased  by  installing  a  larger  electrolyzer.  The 
charging  rate  in  the  second  regime  can  be  increased  by  enabling 
better  cooling  of  the  hydride  tank.  The  transition  point  is  deter¬ 
mined  by  the  electrolyzer  and  hydride  tank  sizes. 

•  System  power  output  during  power  mode  is  governed  at  first 
by  fuel  cell  size  and  later  by  hydrogen  flow  rate.  The  maximum 
power  output  can  be  raised  by  increasing  the  capacity  of  the  fuel 
cell,  whereas  system  runtime  at  a  particular  power  level  can  be 
increased  by  using  a  larger  hydride  tank. 

•  If  integrated  cooling  loops  are  used,  metal  hydride  material 
should  be  carefully  chosen  for  a  particular  application,  not  only  in 
order  to  efficiently  store  hydrogen,  but  also  to  allow  fuel  cell  oper¬ 
ation  at  the  design  temperature  to  enable  operation  at  maximum 
efficiency. 

•  While  major  component  efficiencies  may  improve  for  various 
operating  conditions  during  dynamic  operation,  it  is  important 
to  capture  the  interactions  with  other  system  components  and 
balance  of  plant  items  to  determine  overall  system  performance. 

•  Variable  power  balance  of  plant  components  can  greatly  reduce 
parasitic  losses  for  a  system  with  a  wide  dynamic  operating 
power  range.  Of  course,  variable  power  components  can  increase 
system  cost  and  control  complexity. 
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